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Abstract 

It is shown that after a resummation of leading high-temperature contri- 
butions, a complete and gauge-independent result for the nonabelian Debye 
screening mass at next-to-leading order can be extracted from the static gluon 
propagator. In contrast to previous, incomplete results, the correction to the 
Debye mass is found to be logarithmically sensitive to the nonperturbative 
magnetic mass and positive, in accordance with recent high-statistics results 

from lattice calculations. 
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Over the last few years it has become clear that a complete calculation of perturbative 
corrections to the dispersion laws of quasi-particles in high-temperature QCD requires a 
resummation of the leading-order terms with characteristic scale gT, where g is the coupling 
constant and T the temperature. An improved perturbation theory has been developed 
in particular by Braaten and Pisarski |I| and has been applied in the first place to derive 
damping effects of collective excitations 0, where it turned out to be indispensable in order 
to obtain accurate and gauge independent results for the leading-order term ~ g^T. Most 
recently it has also been employed successfully in a calculation of the next-to-leading order 
term of the QCD plasma frequency Su^^ ~ gu^^ 0. 

All these applications are still within the limits of a perturbative barrier set by the 
presumably incalculable screening of static magnetic fields 0, which causes a breakdown 
of perturbation theory at a certain order. In the case of the gluon self-energy from which 
corrections to dispersion laws are extracted this incalculability is expected at the order 
g^T"^, whereas the above-mentioned results are derived from this quantity evaluated up to 
order g^T"^. However, even at this order occasionally a sensitivity to the behaviour of static 
transverse gluons at momentum scale <^ gT has been found, giving rise to a logarithmic 
enhancement ~ ln(l/5f). 

Besides the dispersion law of propagating collective modes, a quantity of singular impor- 
tance which can be derived from the gauge boson self-energy is the static (chromo)electric 
screening (Debye) mass 0. At leading order IIool^o = 0, k) = rrj-Q ~ g'^T^, which gives rise 
to a pole in the static gluon propagator at k^ = —ml, and, consequently, to exponential 
screening in the potential 
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for gauge group SU(A^) and Nf fermions. In linear response theory, the gradient of $ gives 
the longitudinal electric field generated by a static charge Q. 

Again, one would expect perturbative calculability of the screening mass at the next- 
to-leading order 6 = 5m^ /m^ ~ g. Early attempts (using temporal gauge) ||^,|^ gave 5 = 
—eg, with c a positive constant depending on the particular partial resummation employed, 
whereas in general axial gauge even a different sign was obtained In covariant gauges, 
Toimela found that the time-time component of the gluon self-energy Uoo^ko = 0, k ^ 0) 
is gauge dependent at the order g^T"^, which was confirmed by Nadkarni ||10[ on the basis of 



a dimensionally reduced effective theory. This was taken to mean |]T0[ that the Debye mass 
could not be extracted from the electrostatic propagator. 

However, all these investigations concentrated on the limit k ^ of the static time-time 
component Hoo, guided by the fact that this limit indeed yields the leading-order contribution 
to the Debye mass. But at leading order g'^T'^, Hoo(fco = 0, k) happens to be independent of 
k, whereas the pole of the corrected electrostatic propagator is located at finite |k| ~ gT. 
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Gauge dependence at k ^ 0, away from the physical pole of the propagator, is thus only 
to be expected in a nonabelian gauge theory. Conversely, formal arguments exist that the 
relevant poles of a self-consistently corrected propagator should be gauge-independent also in 
the nonabelian theory And it is just the pole in Eq. (|l|) that determines the exponential 



decay. 

Incidentially, in QED it is equally important to define the electric screening mass in a 
self-consistent manner by the location of the pole rather than the k ^ limit of IIoo as 
almost invariably done in the literature. This in fact modifies the QED Debye mass squared 
at and above the order e^T^ lO . 



According to the resummation program of Ref. a complete calculation of the next- 
to-leading order term in the nonabelian Debye mass should be possible and requires the 
resummation of all hard-thermal- loop contributions. In fact, because only a static quantity 
is to be calculated, this task can be greatly simplified. In the problem of determining next-to- 
leading order corrections to the effective potential in thermal field theories, such a simplified 
scheme has been employed recently by Arnold and Espinosa ||13|. They have found that 



in gauge theories it is algebraically much simpler to resum only static modes. This does 
not touch the completeness of the resummation because, in the imaginary-time formalism, 
nonstatic modes always imply K"^ = — ((27mT)^ + fc^) ^ ^ {qT)"^, so that hard-thermal- 
loop corrections are truly perturbative. Separating the static modes does however give away 
the possibility of a straightforward analytic continuation of any external frequencies. The full 
resummation scheme of Braaten and Pisarski thus is mandatory when external frequencies 
~ gT are to be considered as in dynamical properties of quasi-particles. 

The next-to- leading order contributions to static quantities of relative order 0{g) are 
determined by resummation of one-loop diagrams, and for them it suffices to take into 
account only the static modes (which also precludes fermionic contributions). A further 
major simplification arises in that with all gluon lines being static the hard-thermal-loop 
corrections to the vertices vanish. Only the static gluon propagator is thus needed, which 
reads 



Po=0 
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Po=0 



(3) 



where P = {po, p) and a is the gauge parameter of covariant gauges. 
Evaluating IIoo at relative order 0{g) then yields 



moo(^o = 0, k) = ^-iV / |^{-^ + 1 

;k2 + m2)[3 + 2pk/p2] ^ „ p' + 2pk , 

p (q + m^) P (,q + ^ 




where q = p -|- k and dimensional regularization is understood (yielding no pole terms 
because of the odd integration dimension). Here and in the following the index on m has 
been dropped, the difference in 611 being formally of higher order. 



This result shows that (5noo(fco)k — » 0) is indeed gauge dependent, as found in Ref. |T0 
in a one-loop calculation in the dimensionally reduced effective theory, to which the above 
reasoning has in fact boiled down as far as the 0{g) correction is concerned. It also shows 
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that this gauge dependence is to disappear on an algebraic level when going to the imaginary 
pole = — m^. 

Closer inspection reveals, however, that with — > — there appear "mass-shell" singu- 
larities caused by the massless denominators in Eq. @. A simple introduction of a magnetic 
mass in the denominator in front of the second term of the gluon propagator, Eq. (^, pro- 
vides a physical cut-off, and because the singularities at hand are only logarithmic, the 
coefficient of the ensuing logarithm should be insensitive to the detailed structure of the 



infrared limit of the transverse propagators |14,15|. All the massless denominators figuring 



in the a-independent part of Eq. are indeed associated with the magnetic sector (as 
opposed to pure gauge modes). Assuming mmagn ~ Q^T, the leading contribution to SUqq is 
determined by the logarithmically divergent pieces and reads 



6Iloo{ko,k) 



k2 = 



^ftn^ = 9m^i^^^^l„i + 0(9), (5a) 

giving an unexpectedly large and positive correction to the nonabelian Debye mass. The 
positive sign seems to be a genuine nonabelian effect, for the next-to- leading order correction 



of the Debye mass in e.g. scalar QED is negative [|T^]. Moreover, no logarithmic enhancement 
occurs in the latter. [In spinor QED there are no 0{g) corrections at all due to Pauli 
suppression.] 

The sublogarithmic terms of course do depend on the detailed structure of the infrared 
limit of the transverse propagators, and so cannot be determined completely in the present 
resummed one-loop calculation. However, adopting the hypothesis that this infrared limit 
just amounts to a finite contribution —^]lii(k = 0,k — >• 0) = m'^^g^ ~ g'^T'^ one may go on 
to estimate these sublogarithmic terms from Eq. 

Here one encounters a subtle difficulty with the a-dependent term in Eq. (|^), because 
by approaching the imaginary pole k^ —m?, the explicit factor that apparently ensures 
gauge independence gets cancelled by a linear singularity in the momentum integral. Exactly 
the same phenomenon was encountered in the recalculation of plasmon damping rates in 



general covariant gauges in Refs. , while being absent in homogeneous ( "strict" ) gauges. 
In Ref. I have argued that this behaviour just reflects a singular, gauge dependent 
behaviour of the residue of the propagator rather than an actual gauge dependence of the 
pole determining the dispersion laws. Indeed, introducing an (unphysical) cut-off again 
moves the gauge dependence seemingly afflicting the pole position into the residue, while 
the correction to the pole position becomes independent of this infrared regularization. In 
this way, the sublogarithmic terms are determined |]l9l and give rise to 



5m' = gm'NJ ^ ( In - 1 ) + 0(/). (5b) 

Quite recently, lattice simulations have been performed which permit the extraction of 
corrections to the classical Debye screening length. A rather high precision has been reached 
in simulations of pure SU(2) gauge theory at temperatures up to nearly 8 times the critical 



temperature with the finding of a positive excess in Debye mass squared of 5 = -1-0.30(9) 
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ai P = 3, corresponding to T ^ 7.8Tc and g]^ ~ 1.19. Unfortunately, the coupling is rather 
large so that a quantitative comparison with just the logarithmically enhanced result of 
Eq. ( pal ) is out of question. However, taking the result of Eq. seriously and inserting a 
value for the magnetic mass as suggested by lattice simulations p2|, mmagn ~ 0.1365g^T for 
SU(2), Eq. ( |5ED yields 6 = +0.51, which comes remarkably close considering the largeness 
of g. 

A significantly increased Debye mass has been found previously also in lattice simulations 
of pure SU(3) ||21|, however with larger statistical errors. 

It should be noted that the previous, incomplete results for corrections to the nonabelian 
Debye mass have mostly yielded a negative value for 6, albeit mutually disagreeing in 
numerical magnitude. The complete gauge- independent result presented here differs from 
the former in that all relevant hard-thermal-loop contributions are identified and resummed, 
and also in that the Debye mass is defined through the pole of the electrostatic propagator 
rather than the (gauge-dependent) zero-momentum limit of the gluon self-energy. 

It is the pole that determines the exponential screening in Eq. (|l]), whereas the pre- 
exponential factor therein will generally be gauge dependent. At the next-to-leading order 
considered here, gauge-dependent contributions to the latter indeed arise from the residue 
of the pole of the propagator as well as from a logarithmic branch cut contribution, which 
even gives rise to a (gauge-dependent) modification of the pre-exponential 1/r-behaviour. 
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APPENDIX: 



In this appendix I give the complete contributions to n^i,(A;o = 0,k) at order g^T^ for 
k ^ m gT, in general covariant gauge. The static limit of 11^,^ is transverse, so there are 
exactly two independent structure functions, IIoo and Ua. 

The correction to IIoo as given by Eq. (^) is found to be 



5noo(0,fc) = ^- 

Ztt 
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for k"^ > 0. This can be immediately continued analyticly to k'^ < 0, but for —k^ = m^, 
it has a logarithmic singularity. As argued in the text, the massless denominators in Eq. (^ 
require an infrared cut-off. Let this for simplicity be a mass term A <C m uniformly for 
transverse and gauge modes. In the case of the transverse gluons this can eventually be 
identified with the magnetic screening mass, whereas in the case of the massless gauge 
modes it is just a regulator which drops out from the end results, Eqs. (^). Then for k m 
Eq. ([Al|) gets replaced by 



, , g'^NmT (w? + ^ (m + K 
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In the limit k — > m this gives the gauge-independent result for the Debye mass discussed in 
the text. For general k it remains gauge-dependent, however, so that beyond leading order 
there is no gauge invariant meaning to be attributed to the dielectricity e = 1 — noo(0, k)/k'^ 
as defined from the gluon propagator. 
For k^ 

noo(0,fc)^-i/iVfcT (A3) 

becomes independent of the gauge parameter, but in fact would be different in a background 
covariant gauge p3[. This term would give rise to a modification of the Debye screening in 
a gluon plasma [|^, if such a behaviour persisted for soft momentum k ^ gT. However, the 
complete result of Eq. (AT) reveals that it does not. 



The next-to-leading order correction to the other structure function, njj(0, fc), can be 
derived in a similar manner, and is found to be 



6Yiii{<d,k) = g'NmT 
1 

+ — 
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a + 1)2 + 10 A; 

h 1 -I 

(A4) 
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Again, the magnetic permeability defined by 1/yU = l — ^Uu/k'^ is a gauge-dependent quantity 
beyond leading order. The gauge-dependent terms vanish only at the location of the pole 
of the transverse gluon propagator, which is at k = 0. There the correction term vanishes 
completely, which means that there is no magnetic mass squared of the order g^T"^. The 
magnetic mass must therefore be -C g^^'^T. 

For small k <^ m, Eq. (|A4|) has a linear behaviour with gauge dependent, but positive 
definite coefficient. The transverse propagator therefore has the form l/{k'^ — ck) with 
c ~ g^T. This corresponds to a pole at space-like momentum with k ~ g'^T pSl. However, 



in this regime the perturbatively incalculable contributions ~ g'^T'^ to 11^^ become relevant 
and are expected to remove this pathology by the generation of a magnetic mass term. 
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